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FhaCNDIiG  PATE  H LANK- HOT  ^ILMiD 


I.  INTRODUCTION 


When  the  lateral  surface  of  a long,  unrestrained,  solid  cylinder 
is  exposed  to  a nuclear  thermal  radiation  environment  such  as  that 
indicated  by  Figure  1,  both  a transient,  two  dimensional,  non  axisymmetric 
temperature  field1  and  a transient,  three  dimensional,  thermal  stress 
field  are  produced  within  the  cylinder.  Knowledge  of  the  magnitude  and 
type  of  these  stresses  is  essential  not  only  for  thermal  response  analysis 
but  also  for  the  subsequent  blast  response  analysis.  This  report  describes 
a derivation  of  a set  of  non-dimensional  equations  for  the  stress  field. 
Such  equations  will  provide  a convenient  means  for  calculating  the  values 
of  these  stresses  and  for  performing  a parametric  study  of  the  equation 
variables.  The  results  of  such  a parametric  study  involving  t and  Q, 
two  of  the  three  nuclear  pulse  parameters,  are  also  presented  in  this 
report . 

The  stress  equations  are  derived  within  the  framework  of  classical 
linear  thermoelasticity  theory  using  the  plane  strain  assumption,  and 
the  approach  is  quasi-static  in  which  inertia  effects  and  thermoelastic 
coupling  are  neglected.  Although  the  use  of  the  plane  strain  assumption 
does  not  result  in  a solution  valid  for  all  the  regions  of  the  cylinder, 
it  does  provide  a solution  valid  in  the  regions  away  from  the  ends  of 
the  cylinder. 


II.  TEMPERATURE  EQUATION 

The  equation2  for  the  transient  temperature  field  in  a long,  solid 
cylinder  whose  lateral  surtace  is  heated  by  a nuclear  thermal  environment 
for  which 

1.  the  thermal  properties  of  the  cylinder  are  independent  of 
temperature , 

2.  convection  and  radiation  heat  losses  by  the  cylinder  can  be 
neglected,  and 

3.  the  initial  temperature  field  in  the  cylinder  is  uniform. 


E.F.  Quigley , "Transient  Temperatures  Produced  in  Solid  Cylinders  by 
a Nuclear  Thermal  Pulse",  BRL  Report  No.  1937,  August,  1976,  US  Army 
Ballistic  Research  Laboratory,  Aberdeen  Proving  Ground,  Maryland. 

(AD  HA030703) 
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where  H f(t)  is  the  functional  representation  of  the  curve  in  Figure  1. 
Although  the  temperature  field  is  non-axisymmetric , it  is  readily  seen 
from  (1)  that  the  field  is  symmetric  with  respect  to  $.  Since  the 
temperature  equation  is  needed  for  the  development  of  the  stress  equations, 
it  is  convenient  to  rewrite  (1)  as 
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when  the  initial  temperature  or  the  final  temperature  is  uniform. 

Since  the  initial  temperature  and  the  final  temperature  of  the 
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Y.C.  Fung,  Foundations  of  Solid  Mechanics , Prentic-Hall,  Inc.,  Englewood, 
Cliffs,  New  Jersey,  1965,  p.  C46. 

P.A.  Boley  and  J.H,  Weiner,  Theory  of  Thermal  Stresses,  John  Wiley  & 
Sons,  Inc.,  New  York,  New  York,  1960,  p.  83. 
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Both  the  temperature  and  stress  equations,  (2),  (17),  (18),  (19), 
and  (20)  , can  be  reduced  to  a dimensionless  form  through  the  use  of  the 
following  dimensionless  quantities: 
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IV.  PARAMETRIC  ANALYSIS 


The  values  oT  the  non-diraensional  temperatures  and  stresses  were 
calculated  for  6\=  0.01,  0.1,  and  1.0;  for  t*  = 1.0  and  10.0;  and  for 
9 = 0,  jL,  n , and  £11  . Because  v appears  in  the  right  hand  side  of  (36), 

2 * 

v = 0.33  was  used  ir>  calculating  T Figures  2,3  and  4 are  the  plots 

of  the  temperature  distributions  and  Figures  5 through  14  are  the  plots 
of  the  stress  distributions. 

For  a given  material  of  a given  radius,  the  influence  of  6*  (and 

consequently,  t ) and  Q on  the  temperature  field  can  be  seen  from  Figures 

* 

2,3  and  4.  At  t = 1.0  and  t = 10.0,  the  surface  temperature  for  6=0 

decreases  as  B*  increases  if  Q remains  constant.  For  9 = II  and 

3n  t 

consequent lv , — , the  surface  temperature  increases  as  3*  increases  at 

t*  = 1.0  and  t*  = 10.0.  It  can  be  concluded  from  Figures  3 and  4 that, 
for  S';  1,0,  the  temperature  in  the  cylinder  is  uniform  for  all  values 
of  Q.  For  a given  f*,  it  follows  from  (31)  that  the  temperature  at  any 
point  in  the  cylinder  will  increase  as  Q increases  and  will  decrease  as 
Q decreases. 

The  effect  of  3*  and  Q on  the  thermal  stresses  at  t*  = 1.0  and  jj 
t*  - 10.0  can  be  seen  from  Figures  5 through  14.  Except  for  T*^  on  0=  — 

ami  (See  Figure  8),  an  increase  in  8*  results  in  a decrease  in 

stresses  if  Q remains  constant.  Since  the  lateral  surface  of  the  cylinder 
is  traction  free,  T*^  and  T*-,  are  zero  for  = 1.  Consequently  the 

maximum  value  of  these  stresses  occurs  at  an  interior  point  of  the  cylinde 

as  is  shown  in  Figures  5,  10  and  12.  In  addition  to  influencing  the 

maximum  value  of  these  stresses,  8*  also  effects  the  location  of  the 

maximum  value.  At  t*  = 1.0  and  t*  = 10.0,  an  increase  in  8*  shifts  the 

maximum  value  of  * towards  the  rear  lateral  surface  of  the  cylinder 
rr 

(woe  Figures  5 and  12);  and  an  increase  in  8*  shifts  the  maximum  values 
of  * towards  the  center  of  the  cylinder  (See  Figure  10).  As  in  the 

case  of  the  temperature  distribution,  it  follows  from  (33),  (34),  (35) 
and  (36)  that,  for  a given  6*,  the  stresses  at  any  point  in  the  cylinder 
will  increase  as  Q increases  and  will  decrease  as  Q decreases. 

Using  simple  beam  theory,  this  author  had  shown  that  the  maximum 
axial  stress  is  simply  supported  rectangular  beams  exposed  to  a nuclear 
thermal  environment  did  not  occur  at  t*  = 1.0  hut  at  some  later  t*. 

5 
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By  trial  and  error,  it  was  determined  for  0*  = 0.1  that  maximum  value  of 

T*  , T*  , and  T*  occurred  at  t*  = 1.35  and  that  the  maximum  value  of 
rr  r9’  zz 

T*0  occurred  at  t*  = 1.10.  Figure  15  is  a composite  plot  of  T*f,  T*0, 

T*  and  T*  normalized  to  the  maximum  value  of  T*  , at  the  time  of  their 

zz  r6  zz 

maximum  value.  It  is  seen  from  Figure  15  that  the  maximum  values  of 

T*  and  T*  are  small  compared  to  the  maximum  value  of  T*  . However,  it 
rr  r0  v zz 

is  seen  from  this  figure  that  the  maximum  value  of  T*  is  not  small 

V V 

compared  to  the  maximum  value  of  T*^;  the  latter  being  only  twice  as 

large  as  the  former.  Also,  it  can  be  seen  that  the  axial  stress  is  the 
largest  compression  stress  and  the  largest  tension  stress. 

V.  CONCLUSION 

A set  of  non-dimensional  equations  for  the  quasi-static,  thermal 
stress  field  in  a long,  unrestrained,  isotropic,  homogeneous,  solid 
cylinder  whose  lateral  surface  is  subjected  to  heating  by  a nuclear 
thermal  environment  have  been  derived  under  the  plane  strain  assumption. 
The  results  of  a parametric  analysis  of  the  stress  field  involving  the 
rise  time  of  the  thermal  pulse,  t , and  the  total  energy  of  the  thermal 

environment,  Q,  shows: 

1.  For  a given  Q,  an  increase  in  t results  in  a decrease  in  the 
maximum  value  of  the  stress 

2.  For  a given  Q,  an  increase  in  t shifts  the  location  of  the 
maximum  value  of  the  stress  fo"  those  stresses  which  are  zero 
on  the  lateral  surface  of  the  cylinder,  and 

3.  For  a given  t , the  stresses  at  any  point  in  the  cylinder  will 
increase  as  Q increases  and  will  decrease  as  Q decreases. 
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= Coefficients  defined  by  (29) 

= Dimensionless  variable  defined  by  (31) 

= Coefficients  defined  by  (30) 

= Dimensionless  variable  defined  by  (31) 

= Specific  heat 

= Time  dependent  portion  of  curve  in  Figure  1 
= Radial  coordinate 
= Radius  of  cylinder 

= Dimensionless  variable  defined  by  (31) 

= Time 

= Dummy  integration  variable 

= Rise  time  of  thermal  pulse 

= Dimensionless  variable  defined  by  (31) 

= Coefficients  defined  by  (3)  through  (6) 

= Shear  modulus 

= H f(t) 
o 

= Maximum  irradiance  of  nuclear  thermal  pulse 
= Ordinary  Bessel  function  of  argument  u 
= dJm(u) 


du 


o 


= Temperature  in  cylinder 
= Initial  temperature  in  cylinder 
= Coefficient  of  linear  expansion 
= Dimensionless  variable  defined  by  (31) 

= Angular  coordinate 
= Thermal  conductivity 

= The  n-th  positive  roots  of  ’^(*ranr0)  = ^ 

= I’ossion ’ s ratio 
= Density 
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Dimensionless  variable  defined  by  (31) 

Normal  stresses 
Shear  stress 

Dimensionless  variables  defined  by  (31) 
Temperature  difference 
Solution  of  (11) 

Particular  solution  of  (12) 
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